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ABSTRACT 

We find a new solution of the renormalization group for the Potts model 
with ferromagnetic random valued coupling constants. The solution exhibits 
universality and broken replica symmetry. It is argued that the model reaches 
this universality class if the replica symmetry is broken initially. Otherwise 
the model stays with the replica symmetric renormalization group flow and 
reaches the fixed point which has been considered before. 
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The problem of new critical behavior induced by randomness in spin systems has a 
considerable history. Starting with a classical ip 4 problem, the modified critical behavior 
has later been studied for the two-dimensional Ising and Potts models by various renor- 
malization group techniques, and by numerical simulations. Incomplete list of references 
is provided in [1-12]. Replicas has been used generally to deal with the quenched disor- 
der and replica symmetric solutions have generally been looked for. The first example 
of replica symmetry broken solutions of the renormalization group has been suggested 
in [13], in the context of the p 4 model. In this letter we report on the replica sym- 
metry broken solution for the two-dimensional Potts model with random bonds. The 
model reaches this solution if the replica symmetry is broken initially. In contrast, the 
two-dimensional Ising model turns out to be stable with respect to replica symmetry 
breaking [|J. It reaches always the replica symmetric critical behavior which has been 
studied earlier [5-12]. 

For theoretical study one uses models with a weak disorder, e.g. models with spin 
couplings having small fluctuations around a mean ferromagnetic value. This gives a 
possibility to study the model in continuum, because one reaches the critical point suffi- 
ciently close before the randomness becomes important. For the two-dimensional Potts 
model in particular this allows to use eventually the renormalization group based on the 
conformal theory of the unperturbed model. In this approach the effective theory could 
be described by the Hamiltonian 



where Ho represents, symbolically, the conformal theory of the unperturbed model, while 
the second term with a spatially random mass m(x) coupled to the energy operator 
represents the effective randomness due to spatially inhomogeneous coupling constants 
of spins. Replicating the model and taking the average of the partition function over 
m(x) one gets the effective homogeneous theory with the Hamiltonian: 




(1) 





where g is defined by 




(3) 
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Without loss of generality one could assume a Gaussian distribution for m(x) because 
the terms in the effective Hamiltonian produced by higher moments are irrelevant, in 
the sense of the renormalization group. Also, a spatially uncorrelated distribution of 
coupling constants, and so of m(x), is being assumed. For extra details on the definition 
of this model the reader could consult the papers [8,12]. 

The /3 - function of the renormalization group for the model (2) has been derived, up 
to second order of perturbation theory (third order in g), in [8], with the following result: 

(3(g) = -3eg + 4vr(n - 2)g 2 - 16vr 2 (n - 2)g 3 + 0(g 4 ) (5) 

Here £ is the renormalization group log-scale parameter, and e is related to the central 
charge of the conformal theory for the Potts model at the critical point, on a homogeneous 
lattice. The following parametrization is being used [12]: 

c = 1 - 24a 2 = 1 - 6(a + + a_) 2 (6) 



a± — a ± \J Oq + 1, a + a^ = — 1 (7) 

a 2 + = l + e (8) 

The case of a 2 = |, c = | corresponds to the conformal theory for the Ising model. By 
shifting a 2 + one shifts the central charge c. In particular, for the 3 - components Potts 
model, the parameters of the associated conformal theory are : 

4 „ 6 

c=-, a 2 + = - (9) 

which corresponds to 

e = 4 < 10 > 

The energy operator e of the Potts model corresponds, in the conformal theory classifi- 
cation, to the operator <& lj2 [14]. In general, the dimensions of the conformal operators 
$ n ', n are given by the Kac formula [15]: 

&n',n = «„',„ - 2a «n',n = "j (11) 



For e ~ $1,2 one gets 



1 — n' 1 — n , . 

«n',n = ^ H 2~ a+ ' ^ 



«i,2 = (13) 



A e = 2Ai, 2 = 2(a 2 2 _ 2a «i, 2 ) = -a 2 + - 1 = 1 + -e (14) 

So, for the case of the Ising model, e = 0, A e = 1, the perturbation in (2) is marginal, 
and then one defines the renormalization group for the Potts model in terms of the e - 
expansion technique [8]. 

The renormalization of the operator e(x) has also been found in [8], up to the second 
order, with the following result: 

e(x) -> e ren (x) = Z e (g)e(x) (15) 

^r =lM (6) 

le {g) = 4vr(n - l)g - 8vr 2 (n - l)g 2 + 0(g 3 ) (17) 

Finally, the renormalization of the spin operator has been found in [12], up to the third 
order: 

a(x) -> a ren (x) = Z a (£)a(x) (18) 
d log Z a 



= i«(g) (19) 



lc {g) = -3(n - l)vr 2 e (l + 2 ^Zl|^ri) ) ^ 

+4(n _l) (n _2)7rV + 0(/) (20) 

Here T(z) is the Euler T-function. 

Using results for f3,% and 70-, eqs. (5), (17), (20), where one puts eventually n — 0, 
the following results have been obtained for g c , A £ [8] and A a [12]: 

^-^ + db 2+0 < £3 > < 21 > 

Al = A, - 7e ( 9c ) = A, - |e + + 0( e 3 ) = 1 + g6 2 + 0( e 3 ) (22) 



K = A CT - 7„G/c) = A CT - ^ p( _i^ 2( W^ 3 + 0(e 4 ) (23) 
For the 3-component Potts model, e = — ^, one gets the following numerical values: 

A' £ = 1,02 + (e 3 ) (24) 

A' a = ^- + 0, 00132 + 0(e 4 ) = 0, 13465 + 0(e 4 ) (25) 
15 

In the solution just outlined one assumes that the replica symmetry is not broken 
initially, and then it is preserved by the renormalization group. Physical arguments 
for a more general approach, namely to start with a Hamiltonian in which the replica 
symmetry is lifted, have been suggested in [13]. For the case of the Potts model this 
amounts to replacing the effective Hamiltonian (2) with: 

n „ 

H = ^Hq+ / d 2 x^2 9abSa(x)e b (x) (26) 

a=l J a+b 

where g a b is a Parisi type matrix [16]. We shall assume next this conjecture and give 
the solution of the theory (26), in order to allow, on the basis of the results obtained, to 
verify the conjecture itself. 

By using the techniques of [12] the generalization of the renormalization group equa- 
tions is straightforward. Eqs. (4) and (5) are replaced by 

^ = Pab (27) 
Pab = -^ab + ^(g 2 )ab - 16^ {{g 2 ) aa9ab - { 9ab ) 3 ) (28) 

Here (g 2 ) a b = H c 9ac9cb- The fixed point matrix g ab should satisfy the equation: 

- 3eg ab + 4n(g 2 ) ab - 16n 2 ((g 2 ) aa g ab - (g ab ) 3 ) = (29) 

We always assume the diagonal elements of the coupling matrix g ab to be zero. This 
is because the corresponding terms could always be absorbed into X)" =1 Hq , in eq.(26). 
Then, for the Parisi type matrices, one has in general the following rules [16]: 

9ab -> g(x) (30) 
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(g 2 ) ab - -2g(x) f 1 dyg(y) - f dy(g(x) - g(y)) 2 (31) 

JO JO 

(<7 2 )««-- Cdyg 2 {y) (32) 
Jo 

Here x is the continuous parameter which replaces the matrix indices: g a b ~ g(a — b) ~ 
g(x). Putting r = 3|e| (e is assumed to be negative in its definition (8)), replacing 
g — > ^rg", and using the prescriptions (30) - (32), one gets from (29) the following equation 
for g(x): 

rg{x) - 2gg(x) - f dy(g(x) - g{y)f + g\x) + g{x)g* = (33) 

JO 

Here g = dyg(y), g 2 = dyg 2 (y). Note that the structure of the fixed-point equation 
(33) coincides with the saddle-point equation for the Parisi order parameter function in 
the infinite-range spin-glasses near the phase transition point (the parameter r in eq.(33) 
corresponds to the reduced temperature r = (1— T/T c ) << 1 in the spin-glass model) [16]. 
The solution of this equation is straightforward. Taking a derivative with respect to x 
one gets: 

rg'(x) - 2gg'(x) - 2g'(x) f dy(g(x) - g(y)) 

JO 

+3g'(x)g 2 (x)+g'(x)g 2 = (34) 



So, either 



g\x) = (35) 



or 



r - 2g - 2 f dy(g(x) - g(y)) + 3g 2 (x) + f = (36) 

JO 

Differentiating again one gets 

- 2g'(x)x + 6g'(x)g(x) = (37) 

There are two solutions: 



and 



g(x) = const = gi (38) 



g(x) = ^x (39) 
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Next, we take g{x) in the form: 



< x < Xi 



g{x) = T (40) 
gi, xi < x < 1 



with 



zi = S 9l (41) 
Then, we put back g{x) into the original equation (33). In particular 

g = 9i-\gl g 2 = gl-2gt (42) 

Substituting (40), (42) into (33), either for < x < x\ or for x\ < x < 1 one gets, after 
some simple algebra 

1 1 o 3 3 o , . , 

?i»2 T + f' Xl ~2 T+ 2 ( 3) 

up to the second order in r = 3|e|. This solution can be compared with the replica 
symmetry one, 

g a b ~ g(x) = const = + ^r 2 , < x < 1 (44) 

given by eq.(21) (after a rescaling g ab -> ^^ a6 ). 

Note that in terms of the original "dynamical" eqs.(27),(28) the continuous fixed-point 
solution, eq.(40), is the only one which is (marginally) attractive. This is guaranteed 
by the fact that the eigenvalues spectrum of the corresponding Hessian of the infinite- 
range spin-glass problem is known to be non-positive. On the other hand, all the other 
non-trivial fixed-point solutions which have step-like structure (they correspond to finite 
number of RSB steps in the replica matrix g ab ) are unstable because for any finite number 
of steps there exist positive finite eigenvalues [16]. 

We can now find the dimensions of the operators e and a, for the solution (40). 
Again, with a straightforward generalization, one finds for r y e {g) and 7cr(fi0 the following 
expressions: 

l£ (g) = 4vr^ £ g ab - 8n 2 (g 2 ) aa + 0(g 3 ) (45) 

ab 

iM = 3vr 2 6 (l + 2^r|^y) (g\a + 87r 3 (/) aa + 0(g 4 ) (46) 
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Using the Parisi ansatz g ab — > g(x) with the rules (30)- (32) and, in addition, 




(47) 



one obtains the following expressions: 



i i 




(48) 



o o 




+ Stt 3 / xg 3 (x) + 3g(x) / g 2 (y)dy dx + 0(g 4 ) 




(49) 



Simple analysis of these expressions shows that the modification of A' a for the solution 
(40) for g(x) will be of order e 4 . As we haven't kept the ~ e 4 terms in ^(g) the accuracy 
is not sufficient to use the modified g(x). To e 3 order, A^ remains the same, eqs.(23),(25). 
On the other hand, the accuracy is sufficient for r ) e {g) which depends of terms of order 
~ g and ~ g 2 with coefficients of order ~ 1. Simple calculation with j £ (g) in (45) leads 
to 



which can be compared to the solution of A^ given by eq.(24): A' t = 1, 02 + 0(e 3 ) 

To conclude, we have found an explicit form, at the order e 2 , of replica symmetry 
broken fixed point in the 3-states Potts model with random ferromagnetic bonds. We 
have also calculated one observable quantity, the dimension of the energy operator A", 
which could distinguish this universality class in a numerical experiment [17]. 



K = A £ - le {g) = A £ - \e + 0(e 3 ) = 1 + 0( e 3 ) 



(50) 
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